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Abstract 

The Green's functions are constructed and a one-loop calculation is given for the 
scalar and vector densities and for the equation of state of nuclear matter with a Bose 
condensate of dibaryons. This is the lowest approximation in the loop expansion 
of quantum hadrodynamics, sufficient to account for the presence of dibaryons not 
in the condensate in the heterophase nucleon-dibaryon matter. It leads to a finite 
effective nucleon mass and remains consistent with increasing the density. 
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1 Introduction 



The prospect of observing the long-lived H-particle predicted in 1977 by Jaffe |1[ stim- 
ulated considerable activities in the experimental search for dibaryons 0. Although no 
positive evidence has been observed yet, future experimental studies may elucidate the 
nature of this interesting possibility. Non-strange dibaryons which have a small width due 
to zero coupling to the iViV-channel are promising candidates for experimental searches 
[Q. A method for searching exotic dibaryon resonances in the double proton-proton 
bremsstrahlung reaction is discussed in Ref. [||]. The narrow peaks observed by Bilger 
et al. || and Clement et al. in the pion double charge exchange (DCE) reactions 
on nuclei have been interpreted by Martemyanov and Schepkin J7| as evidence for the 
existence of a narrow d' dibaryon with a mass of 2060 MeV. Recent experiments at TRI- 
UMF (Vancouver) and CELSIUS (Uppsala) seem to support the existence of the d' ||. 
Recently, some indications for a di(1920) dibaryon have also been found ||. 

The existence of dibaryons when settled reliably by experiments should have important 
implications for the properties of nuclear matter. The chemical potential of nucleons 
increases with the density. When it becomes higher than half of the chemical potential 
of dibaryons, the production of dibaryons becomes energetically favourable. Dibaryons 
are Bose particles. In nuclear matter they form a Bose condensate. This interesting 
phenomenon is studied by Baldin et al. \W\ and Chizov et al. [|TT]] by including dibaryon 
interactions through a van der Waals volume correction. A model for nuclear matter 
with an admixture of dibaryons with the short-range nuclear forces approximated by 
a 5-function-like pseudopotential is discussed in Ref. JT^|. The relativistic mean-field 



theory (MFT) with dibaryons was analysed recently by Faessler et al. |13[-||T5[. An 
exactly solvable one-dimensional model for the dibaryon condensation has been found by 
Buchmann et al. [|l6j. Heterophase nucleon-dibaryon matter (HNDM) can be formed 
in interiors of neutron stars |17|]-|]19| and in heavy-ion collisions [15[]. The physics of 
heterophase substances is reviewed in Ref. [p0|| . 

In previous works HIQ"|-|P>|, 113-1231! it was assumed that at zero temperature dibaryons 
in nuclear matter are all in the condensate. In Bogoliubov's model of a dilute interacting 



Bose gas pi] , a fraction of Bose particles is out of the condensate p2| , |23f . This frac- 
tion increases with the density and has an important effect on various thermodynamical 
quantities. 

MFT constitutes the basis of quantum hadrodynamics (QHD) P^]- [p^] . In a loop 
expansion of QHD, it corresponds to the lowest order term (no loops). In MFT all 
dibaryons are in the condensate. Dibaryons that are not in the condensate appear first 
in a one-loop calculation of the scalar and vector densities and of the equation of state 



2 



(EOS). This approximation is identical to the relativistic Hartree approximation (RHA) 
for normal nuclear matter. It is sufficient to account for the existence of non-condensate 
dibaryons. In this paper, we construct Green's functions of HNDM and give a one-loop 
calculation of the scalar and vector densities and the EOS of the binary nucleon-dibaryon 
mixture. 



2 Quantum hadrodynamics with dibaryons 

The dibaryonic extension of QHD is determined by the Lagrangian density 



C = V{id^ -m N - g a a - gjiJfl^ + §(<V) 2 - \m\a 2 - \F 2 V 



(1) 



Here, \I/ is the nucleon field, uj^ and a are fields of the 10- and cr-mesons, = d^oj^ — 
d^ujy is the field strength tensor, and ip is the dibaryon isoscalar-scalar (or isoscalar- 
pseudoscalar) field. Furthermore, and m a are the u- and cr-meson masses and g w , g a , 
hu, hfj are coupling constants of the us- and a-mesons with nucleons (g) and dibaryons 
(h). The w-meson field is described by Stueckelberg's equation |27j. The limit A — > 
recovers the Proca equation. The MFT results [H|-[15] do not depend on the value of A 
in Eq.(®. 

The field operators can be expanded in c-numbers and operator parts: = g^uj c +Cj^ 
a = <r c + a, (p = (p c + (p, and <p* = <p* + (p*. The cr-meson mean field determines the effective 
nucleon and dibaryon masses in the medium: 

m* N = m N + g a a C} 
m* D = m D + h a a c . 

The nucleon scalar and vector densities are defined by the expectation values 

p NS =< M >; (3) 

Pnv =< ^7o^ > • (4) 
The chemical potential equals /ijv = H*n + fc w c an d Hn = ~^*n + Qu^o respectively, for 



nucleons and for nucleon holes in the Dirac sea. Here, pf N = + yW$ + hp is the nucleon 
Fermi energy and kp is the Fermi wave number. 

The dibaryon scalar density is defined as an additional source of the cx-meson mean 
field 

2m* D p DS =< 2(m D + h a a)<p*ip > . (5) 
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The condensate dibaryons give a contribution p c DS = \p c \ 2 to the scalar density. The time 
evolution of the condensate field p c is determined by the chemical potential po- 

<p e {t) = e-^^~. (6) 

To eliminate the time dependence from the condensate parts of the p fields, we pass 
to the /^-representation for dibaryons. This can be done by the substitution p — ► pe % ^ Dt 
and ip* — > p*e~ % ^ Dt . The dibaryon vector density is equal to the timelike component of 
the expectation value of the dibaryon current 

<fi>=<ip*(2ii D -2hdu+i*d) li <p>. (7) 

Here, (^d)^ = (/Ad,0) in the rest frame of the substance. The condensate dibaryons give 
a contribution p c DV = 2p* D p c DS to the vector density with pr> = p* D + h w cu c . 
The self-consistency equation for the effective nucleon mass is: 

m* N = m N - ^r(g a pNS + h a 2m* D p DS ). (8) 

The w-meson mean field is determined from the equation 

rnlujc = gujpNv + h w p DV . 



3 Green's functions of heterophase nucleon-dibaryon 
matter 



The model of Eq.(|T]) is analysed within MFT at zero temperature in Refs. |13|, [14]. To 
go beyond MFT, it is necessary to construct Green's functions of the bosons 

%D AB {x' - x) =< TA{x')B{x) > (9) 

with A, B — a, cj m , (p, p*. For a one-loop calculation, it is sufficient to construct the 
Green's functions in the no-loop approximation. Above the critical density for formation 
of dibaryons, the a- and cj-mesons can be absorbed by the dibaryons in the condensate. 
As a result, the dibaryons leave the condensate and propagate as normal particles. These 
processes occur at tree level, and the Green's functions in the no-loop approximation 
describe the effect of a — uj — tp — p* mixing. 

The Green's functions can be determined self-consistently by solving a system of 
Gorkov-Dyson equations. Multiplying the equation of motion for the a field correspond- 
ing to the Lagrangian density (|I|) by a and taking the time-ordered product, we find the 
average value of the equation over the ground state 
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(-□ - ml) < 7V(1)<7(2) >= 5\1, 2)+g a < T*(1)*(1)<t(2) > 
+2h a < T{m D + h c <j{l))ip*{l)ip{l)a{2) > . 

Taking into account the second-order terms with respect to the operator fields (one-loop 
approximation), one gets 

(_□ _ m l) < Ta(l)a(2) >= 5\l, 2) + 2hlp c DS < Ta(l)a(2) > 
+2m* D K^- s {< T0(l)a(2) > + < T0*(l)a(2) >). 

In the momentum representation, the equation takes the form 



D aa (k) = D™(k) + D™(k)2m* D h a] /p c DS (D^(k) + D* u {k) ) ■ (10) 
The equations for the other Green's functions can be obtained in similar way: 



D™(k) = D™(k)2m* D h^pb s (D^(k) + DP u (k)), (11) 



D^(k) = D™{k)2m* D h^p<h s {D* "(h) + D w (k)), (12) 
D^\k) = D^(k)2m* D h a ^ s (D^(k) + D^\k)), (13) 



D^(k) = D%(k) + D^{k)K^s[{2p* D + k) T Dr(k) + (2p* D - k) T Df (14) 



D^(k) = D^(k)h^p c DS [(2p* D + k) T D^(k) + (2p* D - k) T D* (15) 
Df{k) = D^{k)K^s[{2p* D + k) T D^(k) + (2p* D - k) T D^(k)}, (16) 



£>w (k) — D w (k)+D™ (k)[h^p<h s (2p* D + k) T D^ (k) + 2m* D h a] /p c DS D^ (k)], (17) 
D w \k) = (k)[Ky[fosWD + k) T D^\k) + 2m* D h a Jfa~ s D^(k)\, (18) 



D^'(k) = D^(k)[h U) ^~ s (2p* D - k) T D^"(k) + 2m* D h a ^~D^\k)]. (19) 

Here, 

rr°(k) = t^j, 

k 2 — ml 

f^uiui/i\ _ ~9tiu + kjjtv/m^ _ k^kyjm^ 
ixv v / k 2 -ml k 2 -m 2 ' 

1 



D w (k) 



(k + p* D )*-m% 

are the a- and omeson and the dibaryon MFT propagators. The effective a- and cu-meson 
masses are 

ml = ml + 2hlp c DS , ^ 

m l = m i + 2h ip c Ds- 

The mass of the effective scalar in Stueckelberg's equation is given by m 2 = m 2 w /\. 
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Equations (|T0|) - (|T^) constitute a closed system of equations, which allows us to 
determine selfconsistently the boson propagators including uj — a — tp* — (p mixing. For 
more details see Ref. []TJ] . 



The scalar and vector dibaryon densities are expressed in terms of the dibaryon Green's 
functions 

Dw*(k) = (D^ik)- 1 - ^^(k))/E(k), 

D w (k) = E w (k)/E(k), (21) 

The denominator S(fc) and the self-energy operators are given by the expressions 

S(ife) = {bw*(k)- x - Ew*(it))(^>(fc)- 1 - E^(k)) - E^>*(fc)Sw(fc), 
E^(k) = E^(-k) = {K^-f{2^ D + k)»D"«{k){2f D + k) v 

+ (2m* D h a ^) 2 D-(k), (22) 
E^{k) = = {K^sfi^D + k),D^{k)(2^ D - k) u 

The structure of the dibaryon Green's functions is identical to the structure of the Green's 
functions of fermions in superconductors |23, [28| . 



^From Eqs.(^J) we see that the longitudinal component of the cu-meson propagator 
contributes to the self-energy operators of the dibaryon Green's functions. The model 
is renormalizable for any finite value of A. In the limit A — > the finite renormalizable 
expressions become infinite. Now, we put A = 1 in order to obtain the cu-meson propagator 
used in usual nuclear matter QHD calculations [25, In this case, the QHD results 



for ordinary nuclear matter are reproduced at low densities where there is no dibaryon 
condensate. 

The chemical potential of dibaryons is determined from the relativistic extension of 
the Hugenholtz-Pines relation |29| 



H% - m % = £ w » - S w (0). (23) 

A relation of such a kind is necessary in order to get a pole in the dibaryon Green's 

functions at u = k = and to guarantee the existence of sound in the medium. Eq. (p3|) 
gives 

fa = m* D (24) 

in agreement with MFT [13], ITJ . 
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4 One-loop scalar and vector densities and the EOS 

The one-loop expression for the nucleon scalar density is calculated from Eq.(|3|). After 



renormalization, one gets | 2q , |26 



Pns = Pns + Pns ( 25 ) 

with 

/(^7TO^- |k|) ' (26) 

Pns = -4m 3 N ((m* N /m N ). (27) 

Here, is the statistical factor (7^ = 2 in neutron matter and 7^ = 4 in nuclear matter) 
and 

5 11 
47r 2 C(a;) = x 3 \nx + 1 - x - -(1 - xf + —{I - xf. 

The one-loop expression for the dibaryon scalar density defined by Eq.(|) can be 
written in the form 



1m* D PDS = 2™*dPds + 2m* D < <p*<p > +2h cS J p c DS {< b(p > + < b(p* >). (28) 
We extract the zero order contribution with respect to the dibaryon condensate: 

Pos = ^J-0^D^(k). (29) 
After renormalization, we obtain 

2m *DpDs = ™ 3 D ((m* D /m D ). (30) 
The total dibaryon scalar density has the form 

Pds = Pds + Pds + Pds ( 31 ) 

with 

P " DS = l S 7j|r{^( fc )-^ ( 32 ) 

This integral diverges. The renormalized expression is given by 

1 r+00 

2m * D PDS = ^~2 I dk E kl&E, m* D , p DS ) (33) 

where 

X r ^(k E , m* D , p c DS ) = \ DS (k E , m* D , p c DS ) - X DS (k E , m* D , 0) 

rP d\ DS {k E ,m, D fl) / * \ c d 2 X D s(k E ,m D ,0) 

~PDS dp o Ds \ m D ~ m D)P D S Bm!* D dpP D8 ■ 
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(34) 



We have passed here into the Euclidean space. The density \Ds{kE, m* D , p c DS ) has the 
form 

4(fc|£ w + m* 2 )z 2 + k%{AE u - 3 - k 2 E /m* 2 ) 



\ DS (k E , m * D , P c D5 ) = — , 2 r (35) 



where 



E = kl + mj 
k 2 E + ih 2 J 

z 2 _k 2 E + 8m%R 



R = Pds( 



Am* 2 E w ' 
hi hi 



k\ + rfj^ k% + ml 



The counter terms in the Lagrangian density, responsible for the renormalization of the 
scalar density, are of the form 8C S = (C±a + C2<7 2 )(p*(p. 

The one-loop expression for the dibaryon vector current has the form 



< 3" >= 2(»*d)hPds+ < V*Wd + * d)^ > -2h^p c DS (< uj^ip* > + <u li <p >). (36) 

In the rest frame of the substance, the vector (//£,) M is defined by (p-*d)h = (/•*£>) 0)- 

The nucleon vector density in RHA is the same as in MFT. The dibaryon vector 
density can be represented in the form 



with 



and 



Pdv = Pdv + Pdv + Pdv ( 37 ) 

P ' DV = 1 1 (0 2(/iB + = °> ( 38 ) 

P " DV = l S -0j-A^n + ^)(D^(k)-D^(k)) 

+ ti^y[fo s ) 2 ir*{k){2pr D + u)(V(k) + D^(k))}. (39) 

In this expression, oj is the timelike component of the vector k^. 
The renormalized expression for the value p" DV is given by 

I r+00 

Pdv = J o dk E k 3 E X r ^(k E , m* D , p DV ) (40) 

with 

\™(k E ,m D ,p DV ) = X DV (k E ,m D ,p DV ) - X DV (k E ,m D ,0) (41) 
Here, 

X DV (k E ,m D ,p DV ) = — r . (42) 

2m* D E u ,k%{y/z 2 {z 2 + l) + z 2 ) m D 



8 



The counter term in the Lagrangian density, responsible for the renormalization of the 
vector density, has the form 5C V = C^d^ — ih UJ LO iJ )(p*{d il + ih^u^ip. The structure of the 
counter terms shows that the model is renormalizable to one loop. 

Despite the fact that a dibaryon Bose condensate does not exist in ordinary nuclei, 
dibaryons affect properties of nuclear matter and nuclei through a Casimir effect described 



x) 4 ), so 



by Eq.(0). This effect is analysed in Ref. |(J. When x -»• 1 CO) = 0((1 
the vacuum contributions to the scalar density of nucleons and dibaryons, which have 
opposite signs, are comparable for Ag^/m^ « h^/m^. Dibaryon effects become large 
for h a /{2g a ) m 0.5(4mD/'^Ar) 1 ^ 4 ~ 0.84. Explicit calculation |30| shows that the basic 



properties of nuclear matter at saturation density can be reproduced for h cr /(2g a ) < 0.8, 
and that the results are not very sensitive to the dibaryon mass. The parameter sets of 
the RHA model for the H-particle, and the d'-, and di-dibaryons are shown in Table 1. 
The meson masses m G = 520 MeV and = 783 MeV and the equilibrium Fermi wave 
number kp = 1.3 fm _1 are the same as in RHA without dibaryons [fj^, p6fl . 

The effective nucleon and dibaryon masses, which are obtained as solutions of the 
self-consistency condition (||) are shown in Fig. la. In Fig.l (b-d) the scalar and vector 
densities for nucleons and dibaryons, and the EOS of HNDM are plotted versus the total 
baryon number density pxv for h a /(2g a ) = 0.6 and h^j{2gS) = 0.8. The numerical results 
are given for 77V = 4 (at low density we start from nuclear matter) and for the d' dibaryon 
mass vrtD = 2060 MeV. At densities prv < Pi-, the chemical potential of dibaryons is 
larger than twice the chemical potentials of nucleons, 2p^ < Pd, so we have normal 
nuclear matter. In the interval p\ < pxv < P2, nucleons and dibaryons are in chemical 
equilibrium, i.e. 2p N = p D , and the density of nucleons decreases with increasing total 
baryon number density. In the interval P2 < Ptv < Psi the dibaryon chemical potential 
is inside of the energy gap for nucleons, 2(— m* N + g^oJc) < po < 2{m* N + g w c<j c ), and 
nuclear matter consists of dibaryons only. At densities p% < prv-, antinucleons come into 
chemical equilibrium with dibaryons, the relation 2(—p* N + g^Uc) = po holds true. The 
density of antinucleons increases with the total baryon number density. The numerical 
values for the densities p±, P2, and p% of the nuclear and neutron matter are given in Table 
1. 

The energy density and pressure are calculated from 

£ = J pNdprv and p = J p TV dp N . (43) 
In the interval p 2 < Ptv < P3> one should put 2p N = p D by definition. 
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5 Conclusion 



In an interacting Bose gas, a fraction of bosons is not in the condensate pl|-||23||. Non- 
condensate dibaryons appear in a one-loop calculation of quantum hadrodynamics. This 
approximation corresponds to the relativistic Hartree approximation for normal nuclear 
matter. We have constructed the Green's functions of heterophase nucleon-dibaryon mat- 
ter and have calculated the contribution of non-condensate dibaryons to the scalar and 
vector densities and the equation of state. This contribution increases rapidly with the 
total baryon number density and becomes dominant as compared to the contribution of 
condensate dibaryons. In mean field theory, the effective nucleon mass vanishes with 
increasing density [13|]-[|]1J. In relativistic Hartree approximation, the effective nucleon 



mass is finite and positive. In a loop expansion of quantum hadrodynamics, the relativistic 
Hartree approximation is the lowest approximation sufficient to account for the presence 
of non-condensate dibaryons in heterophase nucleon-dibaryon matter, while providing 
consistent solutions at high densities. 

The present one-loop calculation allows (i) to consider the structure of neutron stars 
with dibaryons in the interiors without restrictions to the total density and (ii) to construct 
the phase diagram for the phase transition of normal nuclear matter to dibaryon matter 
and quark matter without restrictions to the densities and temperatures. 
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Table 1 



Dibaryon Mass (GeV) h a /(2g a ) g c g^ p x j ' p Q p 2 /p p 3 /p ^ N 

3.62 2 

H 2.22 0.6 8.5277 10.4393 11 33 

4.86 4 

2.48 2 

d' 2.06 0.6 8.5438 10.4687 7.9 30 

3.52 4 

1.31 2 

di 1.92 0.6 8.5606 10.4993 4.6 29 

2.23 4 
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Figure captions 



Fig.l. (a) The effective nucleoli and dibaryon masses versus the total baryon number 
density, (b) The scalar density of nucleons of Eq.(|25|): the contribution of (i) nucleons in 
the Fermi sphere (p N g, solid line), (ii) the negative Dirac sea (p' NS , dotted line). The indi- 
vidual contributions to the scalar density of dibaryons according to Eq. (|3T|) : (i) dibaryons 
in the condensate (2m* D p c DS , solid line), (ii) the positive contribution of the zero-point 
dibaryon fluctuations (2m* D p' DS , dotted line), and (iii) the term 2m* D p" DS (dashed line) de- 
scribing the contribution of dibaryons not in the condensate. The densities are measured 
in units of the saturation density po of nuclear matter. Note the different dimensions of 
the scalar nucleon ([fm~ 3 ]) and dibaryon densities ([fm~ 2 ]). (c) The fraction of nucleons 
(pNV I Ptv), an d dibaryons (pdv/Ptv) i n the condensate (solid lines). The fraction of 
dibaryons not in the condensate (p'dv/Ptv) (dashed line), (d) The energy density and 
pressure of heterophase nucleon-dibaryon matter (dashed lines) compared to the RHA 
without dibaryons (solid lines) versus the total baryon number density. At densities 
Ptv < Pi normal nuclear matter (j N = 4) is stable. In the interval p 1 < p TV < p 2 , 
nucleons and dibaryons are in chemical equilibrium. In the interval p 2 < Ptv < Pz, only 
dibaryons are present in nuclear matter. At densities P3 < ptv dibaryons are in chemical 
equilibrium with antinucleons. 
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